clear;clc;
p=[];
n=0;
for c=1:1000
    for k=1:200
        b=c-k;
        a=sqrt(c^2-b^2);
        if(a==floor(a)&&a>0&&b>0)
            if(gcd(a,b)==1&&gcd(b,c)==1&&gcd(a,c)==1)
              fprintf('a= %i,b= %i,c= %i,k= %i\n',a,b,c,k)
              n=n+1;
              p(n)=k;
            end
        end
    end
end
sort(p)